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SUMMARY 

In  order  to  determine  the  applicability  of  the  Zehnder-Mach  inter¬ 
ferometer  to  two-dimensional  cooled-boundary- layer  studies,  an  analyti¬ 
cal  investigation  was  conducted.  It  was  found  that  for  low  wall  to 
free -stream  temperature  ratios,  the  effects  of  light  refraction  may 
cause  considerable  error  in  the  boundary-layer  density  profiles  cal¬ 
culated  from  interferogram  measurements  if  these  effects  are  not  taken 
into  account  and  that  it  is  impossible  to  observe  the  boundary  layer 
very  close  to  the  cooled  wall. 

This  report  therefore  contains  a  discussion  of  the  magnitude  of 
the  refraction  effects  that  may  be  expected  and  acts  as  a  guide  in  the 
design  of  two-dimensional  test  sections  for  the  study  of  cooled  boundary 
layers  by  the  interferometric  method. 

Equations  are  derived  and  used  to  compute  the  trajectories  of  light 
rays  through  particular  boundary  layers  as  well  as  the  approximate  cor¬ 
rections  that  must  be  applied  to  the  boundary-layer  density  values  to 
take  into  account  the  effects  of  refraction.  The  computations  were  made 
for  a  range  of  conditions  that  is  mainly  of  interest  in  turbine-blade 
cooling  studies.  Laminar  and  turbulent  boundary  layers  that  have  a 
free-stream  Mach  number  of  1.0,  no  pressure  gradients  in  the  flow 
direction,  and  wall  to  free-stream  temperature  ratios'  from  approxi¬ 
mately  0.33  to  0.9  were  considered. 


INTRODUCTION 

In  order  to  predict  the  heat-transfer  to  cooled  surfaces  in  hot- 
gas  streams,  local  heat-transfer  coefficients  must  be  estimated. 

Local  heat-transfer  coefficients  cab.  be  obtained  from  interferograms  of 
the  boundary  layer  about  two-dimensional  bodies  in  a  free -convection 
flow  by  determining  the  temperature  gradient  at  the  surface  of  these 
bodies  (reference  l) .  For  forced  convection  of  hot  gases,  the  boundary 
layers  are  usually  very  thin  and  the  temperature  gradients  through  the 
boundary  layer  are  very  large;  light  rays  passing  through  the  boundary 
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layer  -will  therefore  he  refracted  considerably.  Heat  transfer  is  deter¬ 
mined  by  the  region  in  the  boundary  layer  that  is  very  close  to  the 
surface.  When  the  boundary  layer  is  cooled,  light  is  refracted  toward 
the  surface  so  that  the  part  of  the  boundary  layer  nearest  the  surface 
cannot  be  observed.  Even  though  it  is  expected  that  temperature  gradi¬ 
ents  at  the' wall  cannot  be  obtained  in  this  case,  the  prediction  of 
heat-transfer  coefficients  can  be  greatly  facilitated  by  a  knowledge  of 
flow  conditions,  such  as,  how  and  when  transition  to  turbulence  and 
separation  occur  in  the  boundary  layer. 

Corrections  for  the  effects  of  the  refraction  that  are  necessary 
to  obtain  the  true  density  or  temperature  profile  of  a  two-dimensional 
boundary  layer  from  interferograms  are  di&eussed  in  reference  2.  In  a 
talk  on  the  accuracy  of  gas -flow  interferometry  given  at  the  March  5, 
1948  meeting  of  the  Optical  Society  of  America,  F.  Joachim  Weyl  dis¬ 
cussed  the  effects  of  light  refraction  on  the  design  of  two-dimensional 
supersonic  wind  tunnels.  However,  the  case  in  which  light  rays  are 
bent  into  the  surface  was  not  discussed  in  either  of  these  references. 

It  is  desired  to  determine  the  flow  conditions  and  test-section  design 
for  which  the  effects  of  refraction  must  be  considered  and  to  determine 
what  can  be  done  to  minimize  these  effects. 

The  investigation  reported  herein,  which  was  conducted  at  the  NACA 
Lewis  laboratory,  will  aid  in  the  estimation  of  the  order  of  magnitude 
of  the  refraction  effects  that  can  be  expected  in  an  interferometric 
investigation  of  two-dimensional  boundary  layers  similar  to  those 
encountered  in  the  flow  over  cooled  turbine  blades.  Numerical  results 
are  obtained  for  flow  in  laminar  and  turbulent  boundary  layers  on  a 
flat  plate  at  a  free -stream  Mach  number  of  1.0  and  at  wall  to  free- 
stream  temperature  ratios  from  approximately  0.33  to  0.9. 

METHOD  OF  ANALYSIS 

The  interferometer  is  best  suited  to  two-dimensional  studies > 
therefore,  the  boundary  layer  existing  on  a  cooled  model  contained  in 
a  rectangular  channel  will  be  considered.  It  is  assumed  that  the  model 
spans  the  distance  between  two  glass  walls  and  has  its  stagnation  point 
at  some  point  downstream  of  the  entrance  to  the  channel. 

It  is  also  assumed  that  various  conditions  of  pressure,  gas  veloc¬ 
ity,  and  temperature  can  be  attained  at  any  point  in  the  channel.  The 
boundary  layer  to  be  investigated  is  initially  laminar  and  then  transi¬ 
tion  to  turbulence  takes  place.  According  to  reference  3,  transition 
can  transport  itself  laterally  across  the  channel  by  the  mechanism  of 
transverse  contamination.  In  such  a  case,  the  transition  region  is 
V-shaped.  The  glass -wall  boundary  layer  and  the  boundary  layer  on  the 
model  surface  intersect  at  the  corners  of  the  channel  as  shown  in 
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figure  1,  which  is  a  .schematic  diagram  describing  the  path  a  light  ray 
would  take  when  passing  through  a  boundary  layer  formed  in  a  "two- 
dimensional"  test  section  during  forced-convection  flow  over  a  cooled 
surface.  The  free -convection  boundary  layers  on  the  outside  of  the 
channel  due  to  the  temperature  differences  between  the  glass  wall  and 
outside  air,  the  nonuniform  temperature  distribution  in  the  glass  walls 
due  to  heat  flow  to  the  cooled  model,  and  the  inside -wall  boundary 
layers  as  well  as  the  boundary  layer  under  investigation  affect  the 
trajectory  of  a  light  ray  passing  through  the  channel.  In  this  investi¬ 
gation,  it  is  assumed  that  the  outside  free  convection  has  been  effec¬ 
tively  eliminated  and  that  the  density  gradients  in  the  glass  walls  and 
glass -wall  boundary  layers  are  parallel  to  the  ray  entering  the  glass 
wall,  A  diagram  of  a  test-section  design  that  is  expected  to  allow  the 
previous  assumptions  to  be  approximately  satisfied  is  presented  in  fig¬ 
ure  2.  The  cooled  body  is  removed  from  the  glass  wall  in  order  to 
reduce  the  local  thermal  variation  of  the  index  of  refraction  in  the 
glass  wall  and  to  prevent  the  intersection  of  the  boundary  layers  at  the 
corners  of  the  channel.  The  effect  of  the  gap  size  between  the  model 
and  the  glass  on  the  apparent  thickening  of  the  boundary  layer  has  not 
been  investigated,  but  it  is  known  that  removing  the  cooled  body  from 
the  window  considerably  reduces  the  thermal  distortion  of  the  glass. 

The  evacuated  chambers  on  the  sides  of  the  channel  allow  thin  glass  to 
be  used  in  the  channel  wall  when  the  air  flow  in  the  channel  is  at  low 
pressure  and  eliminate  much  of  the  free  convection  on  the  outside  of 
the  channel.  The  evacuated  chambers  require  that  glass  be  inserted  in 
the  reference  path  of  the  interferometer  to  compensate  for  the  thick 
glass  on  the  chambers.  1 

The  two  effects  of  refraction  that  will  cause  difficulties  in 
evaluating  interf erograms  are:  (l)  The  lower  light  rays  are  bent  into 
the  surface  and  therefore  will  not  appear  at  all  on  an  interf erogramj 
(2)  Because  a  light  ray  passes  through  regions  of  varying  density,  a 
fringe  shift  appearing  on  an  interferometer  screen  indicates  a  value  of 
the  density  integrated  on  a  curved  light -ray  path  through  the  boundary 
layer . 

The  first  effect  is  investigated  by  obtaining  an  expression  for  the 
trajectories  of  light  rays  through  various  laminar  and  turbulent  bound¬ 
ary  layers.  This  expression  is  presented  in  a  dimensionless  form  that 
avoids  numerical  calculations  involving  various  Reynolds  numbers,  free- 
stream  densities,  light  sources  of  different  wavelengths,  and  distances 
back  from  the  stagnation  point.  The  form  of  the  expression  is  such  # 
that  numerical  representations  of  boundary- layer  temperature  profiles 
may  be  readily  substituted. 

4 

All  symbols  axe  defined  in  appendix  A;  the  equation  of  the  trajec¬ 
tory  of  a  ray  is  given  in  appendix  B  as: 
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where 


(\r\ 


(1) 


and 


The  height  parameter  tj  is  essentially  the  ratio  of  the  distance  from 
the  surface  y  to  the  houndary-layer  thickness  5.  In  the  expression 
for  the  parameter  tj  it  is  assumed  that  the  houndary-layer  thickness 
is  proportional  to  some  power  of  x.  This  kind  of  behavior  is  quite 
often  found  in  laminar  and  turbulent  boundary  layers. 

The  second  effect  can  he  estimated  by  considering  that  a  light  ray 
entering  the  channel  at  position  y-j_  (see  fig.  3)  appears  to  traverse 
the  channel  at  a  height  ya  and  indicates  a  density  pa  on  an  inter- 
ferogram.  The  second  effect  can  therefore  he  described  by  the  ratios 
pa/pi  arid  ya/yi*  The  ray  that  passes  through  a  point  at  a  height  y 

in  the  test  section,  reaches  an  image  point  on  the  interferometer 
screen.  The  point  at  a  height  ya,  or  position  A  in  the  object  plane, 
corresponds  to  this  image  point.  In  appendix  B  it  is  shown  that 

«•■  I2e2)  <*> 

where 

S  .  _ (I 

K(kpg)  Be”  (L'/x) 

'  '  '  l 

The  equation  used  to  evaluate  int erf erograms  where  the  refraction 
effects  are  small  is  .  v 

sx  =  L  (na  -  nQ)  (3) 

By  substituting  - 

n  s*  kp  +  1  (4=) 
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auantitvlndeX  °f  rffr®;ctlon  T  ln  equation  (3)  an  expression  for  the 
quantity  pa  can  he  found.  The  value  of  the  fringe  shift  S  at  the 

image  point  which  corresponds  to  ya  in  the  object  plane  is  to  be  sub 

stituted  in  the  equation  in  order  to  obtain  pa.  In  appendix' B,  the 

that°takes  ^tnn§th  S+  +JS  alS°  obtained  an  integration  procedure 

equatioLeSfSt0SACCa0rTcoiein;draCti0n  °f  **  ^  ray*  When 


1/1-2 


1/Q 


SrLnSiin?13  0n+the  r/Sht  Side  of  the  equation  axe,  respectively, 
the  contributions  to  due  to  the  boundary  layers  on  the  glass 

waHs  inside  the  test  section,  the  refraction  of  a  ray  through  the  ■ 
boundary  layer  under  investigation,  and  the  contribution  due  to  the 

SSTfcSJLSV  point  s the  lMide  °f  —t  wind*  L  Wf Lge 

Plane  (reference  2) .  The  three  contributions  were  found  to  be  6 

■  (tz\  2!i8f 

Vpi/f  rf  L  (6) 


0  a/T  dS 


1/1-2 


2a/^  \J  0  r  aJ\ 


* 

mS\oSda£  ValUe  °f  tie  ********  »tlo  In  the  glass 


NUMERICAL  calculations 

is  desired  to  demonstrate  the  applicability  of  the  nrecedinv 
analysis  in  estimating  the  magnitude  of  the  refraction  effects  tha/are 
encountered  In  an  Interferometric  Investigation  of  Cn"  2d 
lent  boundary  layers  similar  to  those  encountered  on  cooled  turbine 
blades  and  in  designing  a  test  section  that  will  be  useful  in^bSinin, 
heat-transfer  information  from  interferograms  of  the  bwSa£  Syer?  ® 

over  all^S^^T  described  in  this  section  are  for  flow 
flat  plate  at  a  free-stream  Mach  number  of  1.0.  However,  it  is 


* 
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possible  to  use-the  method  presented  for  flows  in  which  thereare  pres¬ 
sure  gradients  present  and  that  have  Mach  numbers  °therthan  1.0^  It 
is  not  the  purpose  of  this  report  to  describe  the  effects  of  light 
refraction  for  all  free -stream  conditions  that  may  be  encoun  ere 
rather  to  show  general  trends  that  may  be  expected. 


Temperature  Profiles 

The  theoretical  temperature  profile  over  a  flat  plate  for  a  com¬ 
pletely  laminar  boundary  layer,  as  obtained  from  reference  4  and  from 
additional  information  supplied  to  the  MCA  by  Dr.  Crocco,  is  presented 
in  figure  4(a).  The  ratio  of  the  absolute  static  temperature  in  the 
boundary  layer  to  the  absolute  free-stream  static  temperature  is  • 
plotted  against  the  usual  parameter  r\.  The  property  va  ues  ' 

obtaining  the  boundary-layer  temperature  profiles  are  f ^airandare 
assumed  to  vary  as  a  power  function  of  the  temperature.  The  functi 
used  for  the  viscosity-temperature  and  specific  heat -temperature  rela- 

(  T  \°'75  d  °P - (jL^  >  respectively.  The 

tions  are  £  =  cp,&  \TbJ 

Prandtl  number  is  equal  to  0.725.  For  laminar  boundary  layers,  the 
exponent  of  the  Reynolds  number  is  m  =  0.5  and  K  =  1 
expression  for  tj. 

In  order  to  obtain  a  corresponding  plot  for  the  turbulent  boundary 
Inver  the  similarity  that  exists  at  Pr  =  1  between  the  velocity 
field' and  the  field  of  total  temperature  drop  (reference  5)  is  used  to 
obtain  an  expression  for  r  in  terms  of  T). 


T’  -  T 


w 


T’  -  T 
&  w 


(9) 


By  substituting 


'  =  t(i  +  ^  M2) 


(10) 


and 


M  =  Mg 


u  (T 

U  ^  r 


(ID 


for  the  total  temperature  and  Mach  number  in  equation  (9)  it  can 

shown  that 

r  =  ^  %2  §  i1  "  S)  +  rw  l1  “  §)  +  U 


he 


(12) 
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In  the  turbulent  region  of  the  turbulent  boundary  layer  the  velocity 
profile  is  approximated  by  a  power  law 


u  /jf/H 
U  “  \b) 


(13) 


Measured  velocity  profiles  for  low  subsonic  flow  show  that  the  value 
of  N  can  be  taken  as  seven.  In  reference  6,  it  is  shown  that  the 


boundary  layer  thickness  is 

•n 

v  Re 

Then  rj  =  £  =  K -  y  where 

o  X 


8  -  °'376  X 

S  -  “1/5  x 
Re  ' 

K  =  2.66, 


if  the  l/7  power  law  is  used, 
and  m  =  0.2. 


In  the  laminar  sublayer,  it  is  assumed  that 


T  -  T. 


w 


T  - 

s 


w 


(14) 


and  in  reference  6 


^S 

5 


191 


Re 


0.7 


(15) 


for  a  turbulent  boundary  layer  following  the  l/7  power  law.  From 
equations  (12),  (14),  and  (15),  it  can  be  shown  that 


r  =  constant  r|  +  rw  (16) 

in  the  laminar  sublayer.  The  temperature  ratio  for  the  laminar  sub-, 
layer  and  the  turbulent  region  are  plotted  in  figure  4(b)  for  a  Mach 
number  of  1.0.  For  the  laminar  sublayer  the  Reynolds  number  was  chosen 
at  5X10^  in  order  to  obtain  the  constant  in  equation  (16).  The  turbu¬ 
lent  portion  of  the  profiles  described  by  equation  (12)  is  independent 
of  Reynolds  number  as  long  as  the  l/7  power  law  is  valid. 


Light-Ray  Trajectories 

The  possibility  of  calculating  trajectories  of  light  rays  through 
the  boundary  layer  is  offered  by  equation  (l).  In  figure  5,  the  light 
paths  through  the  three  laminar  and  the  three  turbulent  boundary  layers 
shown  in  figures  4  and  5  are  described.  The  right  side  of  equation  (l) 
*  was  evaluated  by  means  of  Simpson^  rule.  The  slope  of  the  curwes  in 

figures  6  and  7  depend  only  upon  the  parameters  T)  and  rv  as  shown 
.  in  equation  (l).  Therefore,  for  any  value  of  r\  all  curves  for  the 
same  parameter  rv  have  the  same  slope.  This  fact  is  helpful 
in  interpolating  .values  between  the  curves. 


I 
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Determination  of  Ygjvi 

From  equation  (2),  is  against  t2  for  various 

values  of  rv  and  t|-]_  in  figure  6.  It  is  assumed  that  the  object 

plane  corresponds  to  the  center  plane  of  the  channel,  that  is, 

€ «  =  1  a  C. 


Determination  of  p  /p, 

#  ^ 

An  estimate  of  the  quantity  p  /p,  for  the  prescribed  boundary 
layers  may  be  made  from  figures  7  to  9.  In  figure  7,  (Pa/Pp)p  is 
plotted  against  Sf/L  for  various  values  of  the  quotient  r-Jxf. 

Plots  of  (Pa/pi)i-2  against  |2  for  various  wall-to- stream  tempera- 
ture  ratios  and  values  of  are  presented  in  figure  8.  The  calcula¬ 

tions  for  these  plots  were  made  by  assuming,  in  equation  (7)  ,  that  the 
glass -wall  boundary  layer  is  negligible,  that  is,  Sf/L  =  0.  For 
values  of  6f/L  >0  the  ratio,  (Pa/Pl)  1-2  should  be  multiplied  by 

^1-2  -~j .  The  values  of  (pa/pi)i_2>  8X6  limited  because  the  rays  are 

bent  into  the  surface. 

In  figure  9,  (Pg/Pl) ft  .  is  plotted  against  |2  for  various  tern- 

perature  ratios  and  values  of  T]^.  The  ob  ject  plane  was  assumed  to 
correspond  to  the  plane  containing  the  center  plane  of  the  test  section. 


A  Numerical  Example 


It  is  desired  to  design,  a  channel  that  can  he  used  to  obtain 
local  heat -transfer  coefficients  over  a  flat  surface  by  means  of  the 
Zehnder -Mach -interferometer  for  a  given  air  temperature,  wall-to- 
stream  temperature  ratio,  Reynolds  number,  and  Mach  number  of  the  order 
of  magnitude  of  1.0,  The  design  procedure  consists  of  assuming  various 
values  of  the  pressure  and  temperature  of  the  air  in  the  test  section 
and  then  determining  if  reasonable  values  of  the  fringe  shift  across 

x  (Pa/Pl)f 

the  boundary  layer,  the  ratio  7*,  and  the  ratio  - / *■—  are 

L  V°1 

obtained.  The  effects  of  refraction  that  can  be  expected  for  a  par¬ 
ticular  design  must  be  estimated  before  these  calculations  may  be  made. 


In  order  to  obtain  the  local  heat -transfer  coefficient  directly 
from,  the  fringe  shifts  on  an  interf erogram,  it  is  necessary  to  calculate 
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the  boundary- layer  temperature  gradient  at  the  wall.  However,  all  light 
rays  are  bent  toward  the  'surface  and  it  will  not  be  possible  to  view  the 
boundary  layer  near  the  wall.  For  this  reason,  it  is  assumed  sufficient 
to  obtain  the  temperature  profile  up  to  the  point  at  which  the  tempera¬ 
ture  gradient  is  0.9  of  the  gradient  at  the  wall  in  the  laminar  boundary 
layer  and  up  to  a  value  of  y/&  =  0.01  in  the  turbulent  boundary  layer. 
The  profile  shapes  shown  in  figure  4  were  used  to  determine  the  corre¬ 
sponding  value  of  tj  in  the  present  evaluation.  In  addition,  it  has 
been  found  that  it  is  very  difficult  to  obtain  the  temperature  profile 
for  points  that  are  very  close  to  the  wall.  The  boundary- layer  thick¬ 
ness  is  therefore  made  large  enough  so  that  the  portion  of  the  boundary 
lg.yer  which  vanishes  is  at  least  0.004  inch. 

Because  it  is  desired  to  keep  the  glass -wall  boundary  layer  thin 
it  would  be  helpful  to  so  design  the  rectangular  channel  that  the 
length  of  the  wall  containing  the  glass  is  as  short  as  possible.  This 
design  is  accomplished  by  using  one  wall  of  the  channel  as  the  test 
surface,  which  is  considered  to  be  cooled  over  its  entire  length. 

Air  at  70°  F  is  assumed  to  flow  through  the  channel  at  a  maximum 
Reynolds  number  of  5X10^.  The  pressure  in  the  channel  Is  maintained 
constant  at  about  26  inches  of  mercury  vacuum.  A  test-section  design 
is  obtained  for  three  wall-to-stream  temperature  ratios  for  both  laminar 
and  turbulent  boundary  layers.  The  light  source  in  the  interferometer 
is  a  filtered  mercury  arc  lamp  giving  off  only  the  green  line  which  has 
a  wave  length  of  5461  angstrom  units. 

The  results  of  the  calculations  are  presented  in  table  I.  In 
column  1  of  this  table,  boundary  layers  having  three  different  wall-to- 
stream  temperature  ratios  each  are  selected.  In  column  2,  the  minimum 
values  for  the  distance  back  from  the  stagnation  point  are  obtained  by 
substituting  the  lowest  value  of  tj-,  that  is  desired  into  the  relation 
Rem 

tj  =  K - y  where  y  is  set  equal  to  0.004  inch  and  the  Reynolds 

X  5 

number  is  the  maximum  value,  5X10  .  The  maximum  values  for  the  span, 
column  3,  are  obtained  from:  the  value  of  from  figures  5(a)  and 

5(b),  which  are  plots  of  the  trajectories  of  light  rays  through  the 
boundary  layer;  the  maximum  Reynolds  number;  a  value  of  kp^  corre¬ 
sponding  to  the  green  line  from  a  mercury  arc  lamp;  a  pressure  of 
26  inches  of  mercury  vacuum;  an  air  temperature  of  70°  F;  and  the 
value  x,  in  column  1.  The  ratio  z/l1  is  equal  to  1.0  and  the  maxi¬ 
mum  velocity,  in  column  4,  is  obtained  at  the  maximum  Reynolds  number. 

In  column  5,  values  of  €  at  ya/y-^  =  1  are  presented;  C*  is 

related  to  €  by  the  relation 

€L/2  -  S„ 
e'  =  r 
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where  Sf  Is  the  glass-wall  h oundary- layer  thickness  found  hy  assuming 
the  wall  boundary  layer  to  be  nearly  isothermal  and  to  have  Its  stagna¬ 
tion  point  at  the  same  point  as  that  for  the  cooled  boundary  layer. 

The  values  of  p  /pn  are  found  from  equations  (5)  to  (8)  where  rf  is 

assumed  equal  to  1  because  the  glass-wall  boundary  layers  are  nearly 
isothermal  and  T]_  is  the  value  of  r  corresponding  to  Column  6 

gives  the  differences  between  the  density  at  position  one  and  the 
apparent  density  as  compared  with  the  density  at  position  one.  In 
column  7,  the  fringe  shift  that  will  appear  on  an  interferogram  between 
the  free  stream  and  the  wall  is  given.  The  equation  for  the  fringe 
shift  from  the  free  stream  to  any  point  in  the  boundary  layer  is  derived 
from  equations  (3)  and  (4)  and  is 

S»-6  -  x  <  V  <17) 

This  equation  is  plotted  in  figure  10  for  the  conditions  of  this  problem 
a.nfl  for  a  span  of  1  inch.  A  positive  sign  for  the  fringe  shift  indi¬ 
cates  that  pa  is  greater  than  pg.  This  fact  is  discussed  in  the 
next  section  of  this  report .  The  ratios  of  the  minimum  distance  x  to 
the  maximum  span  is  given  in  column  8.  In  column  9  ,  the  ratios  of  the 
portions  of  Pa/p^  that  result  from  the  glass-wall  boundary  layers  to 

the  total  9qJp\  are  listed. 


DISCUSSION  OF  RESULTS 

The  principal  factor,  other  than  the  flow  conditions,  governing  a 
test  section  design  is  the  estimate  of  how  close  to  the  wall  the  bound¬ 
ary  layer  need  be  studied;  that  is,  a  minimum  value  of  y/b  must  be 
chosen.  In  some  cooling  studies,  it  may  only  be  necessary  to  deter¬ 
mine  whether  a  particular  portion  of  a  boundary  layer  is  laminar  or 
turbulent  or  where  separation  has  taken  place.  These  types  of  investi¬ 
gation  should  be  comparatively  simple  with  the  interferometer.  For 
some  investigations,  however,  it  is  necessary  to  obtain  the  temperature 
gradient  at  the  wall.  Because  the  temperature  ratio  is  less  than  1.0 
and  all  rays  are  bent  toward  the  surface,  there  is  a  definite  limit  as 
to  how  close  to  the  surface  the  boundary  layer  may  be  investigated. 

The  restriction  on  the  minimum  value  of  y/b,  for  which  it  is  desired 
to  obtain  the  temperature  in  the  boundary  layer,  also  restricts  the 
boundary-layer  thickness  6  to  some  minimum  value  because  there  is  a 
-minimum  value  of  y  for  which  the  profile  can  be  investigated. 

In  general,  there  are  four  effects  which  limit  the  minimum  y 
value. 
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(1)  If  at  best,  it  is  possible  to  estimate  to  l/lO  of  a  fringe  in 
determining  the  fringe  shift,  the  number  of  fringe  shifts  across  the 
boundary  layer  limits  the  minimum  y  value  that  can  be  investigated. 

(2)  Even  with  no  flow  over  a  surface  a  slight  fringe  shift  is  often 
obtained  close  to  the  surface.  This  shift  may  be  due  to  reflection  of 
light  from  the  surface.  It  is  conceivable  that  this  reflection  limits 
the  accuracy  with  which  fringe  shifts  can  be  estimated  close  to  a  sur¬ 
face  as  a  result  of  interference  between  the  survey  rays  and  the 
reflected  rays. 

(3)  The  effect  of  diffraction  of  light  from  the  edge  of  the  model 
results  in  diffraction  fringes  parallel  to  the  model  surface  which 
cause  difficulty  in  determining  the  interference  fringe  shift. 

(4)  There  is  also  some  difficulty  in  obtaining  the  exact  location 
of  the  model  surface  on  an  interferogram  due  to  the  previously  men¬ 
tioned  diffraction  effect)  the  intersection  of  fringes  with  the  surface 
therefore  cause  some  error  in  "the  measured  positions.  All  these  effects 
contribute  to  the  great  difficulty  encountered  in  obtaining  the  tempera¬ 
ture  profile  in  the  boundary  layer  very  near  to  a  model  surface. 

The  ratio  of  the  distance  from  the  stagnation  point  x  to  the 
span  of  a  channel  L  is  an  indication  of  the  feasibility  of  a  certain 
test-section  design.  For  large  values  of  the  ratio  x/L,  it  can  be 
expected  that  end  effects  such  as  the  glass -wall  boundary  layer  and 
transverse  contamination,  in  the  case  of  laminar  boundary  layers,  will 
be  larger  and,  accordingly,  will  influence  the  interferogram  fringe 
shift  to  a  considerable  extent.  It  is  difficult  to  evaluate  the 
average  density  in  these  end-effect  regions,  therefore  it  is  reasonable 
to  believe  that  it  would  be  difficult  to  evaluate  interferograms  for 
large  values  of  x/L.  This  observation  is  illustrated  in  columns  8  and 
9  of  table  I.  The  contribution  to  p  /p  due  to  the  glass -wall  bound- 

3/  J-  p 

ary  layers  increases  as  the  ratio  x/l  increases. 

The  value  of  e,  for  ya/y-]_  =  1,  varies  but  little  from  0.8  for 
the  range  of  conditions  investigated  (table  I,  column  5). 

i 

An  evaluation  procedure  based  on  the  assumption  that  the  refrac¬ 
tion  is  small  would  cause  considerable  error  in  the  profiles  calculated 
from  fringe  shifts  as  evidenced  by  column  6  of  table  I.  The  negative 
fringe  shift  in  column  7  indicates  that  it  would  be  possible  to  obtain 
a  wall  temperature  greater  than  the  free-stream  temperature  if  a  cor¬ 
rection  were  not  made  for  refraction.  This  result  may  be  in  error 
because  terms  of  higher  degree  than  the  second  in  6  were  neglected. 

The  results  shown  in  column  7,  however,  still  show  that  the  effects  of 
refraction  increase  with  decreasing  temperature  ratio,. 
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For  given  values  of  “  |  K2  (kp&)  Re  Reynolds  number, 

and  gas  temperature}  tlie  value  of  x/l  can  decrease  for  decreasing 
pressure.  F uriiiermore ,  for  a  constant  Reynolds  number  and  gas  tempera¬ 
ture,  a  very  low  pressure  is  associated  with.  a  high  velocity*  These 
conditions  of  low  pressure,  high  velocity,  and  small  x/L  are  encoun¬ 
tered  in  hypersonic  air -flow  channels .  However,  with  the  very  low 
pressures  encountered  in  high  Mach  number  flow  it  can  be  expected  that 
smai  1  fringe  shifts  will  occur  across  the  boundary  layer. 

The  results  of  a  preliminary  investigation  for  the  hypersonic  range 
indicate  that  the  interferometric  method  will  not  he  useful  as  a  means 
I  of  investigating  laminar  boundary  layers  for  temperaure  ratios  near  one 
because  of  the  very  small  fringe  shift  across  the  boundary  layer  but 
will  be  useful  for  laminar  layers  with  low  wall-to-stream  temperature 
ratios.  In  the  case  of  turbulent  boundary  layers,  investigations  in  the 
range  of  temperature  ratios  near  one  are  also  ruled  out  because  of 
insufficient  fringe  shift  across  the  boundary  layer  and  the  low  wall- 
to-stream  temperature  ratio  is  ruled  out  because  the  effects  of  refrac¬ 
tion  are  expected  to  make  an  accurate  evaluation  of  the  interf erogram 
impossible. 


to 

£ 


CONCHISIOHS  AHD  KECOMMEHDATIONS 

A  method  has  been  presented  for  estimating  the  effects  of  light 
refraction  in  an  interferometric  investigation  of  two-dimensional 
boundary  layers. 

From  plots  of  the  trajectories  of  light  rays  through  both  laminar 
and  turbulent  boundary  layers,  from  calculations  of  the  error  involved 
in  not  taking  into  account  refraction,  and  from  the  results  of  a  numer¬ 
ical  example}  it  was  found  that  the  effects  of  refraction  become  more 
prc^LOunced  as  the  wall— to— stream  temperature  ratio  is  decreased. 

The  interferometric  method  is  believed  to  be  limited  as  a  method  of 
obtaining  local  heat-transfer  coefficients  for  forced-convection  flow  at 
low  wall-to-stream  temperature  ratios.  The  principal  cause  of  the 
limitation  on  the  method  is  the  large  temperature  gradient  which  causes 
rays  to  be  bent  into  the  surface  of  the  model  and  prevents  the  observa¬ 
tion  of  that  portion  of  the  boundary  layer  very  close  to  the  wall.  In 
addition,  surface  reflection,  diffraction,  and  uncertainty  in  locating 
|  the  model  surface  on  an  interf  erogram  limits  measurements  on  inter - 
ferograms  very  near  the  wall. 
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It  is  believed  that  the  interferometer  would  he  useful  for  investi¬ 
gating  laminar  boundary  layers  in  hypersonic  channels  with  low  wall-to- 
stream  temperature  ratios  as  well  as  low  velocity  flows  with  wall -to - 
stream  temperature  ratios  near  one. 

Removal  of  the  cooled  model  from  the  glass  wall  would  reduce  or 
eliminate  the  apparent  thickening  on  interferograms  owing  to  the 
intersection  of  the  model  boundary  layer  and  the  glass-wall  boundary 
layers  at  the  concave  corners  of  a  channel  and  reduce  or  eliminate  the 
thermal  distortion  of  the  glass  walls  that  would  exist  if  the  model 
contacted  the  glass. 

The  glass  wall  should  be  made  as  thin  as  possible  so  that  the 
refraction  effects  due  to  any  thermal  distortion  of  the  glass  may  be 
neglected. 

In  order  to  allow  thin  glass  windows  to  be  used  with  a  low  pres¬ 
sure  in  the  channel  it  would  probably  be  necessary  to  provide  evacuated 
chambers  on  the  outside  of  the  channel.  The  evacuated  chambers  would 
have  the  additional  advantage  of  reducing  free  convection  on  the  glass 
test-section  walls. 


Lewis  Plight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  May  17,  1951. 
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APPENDIX  A 


SYMBOLS 


The  following  symbols  are  used  in  this  report: 

Cp  specific  heat  at  constant  pressure,  Btu/(lb)(°F) 

cv  specific  heat  at  constant  volume,  Btu/(lb)(°F) 


L  T 

11  €2  “  3>  ft 


K  constant  in  the  expression  r\  =  K  y,  dimensionless 

k  Dale -Gladstone  constant,  cu  ft/ slug 

L  span,  inside  test  section,  ft 

L'  L  -  25f,  ft 

M  Mach  number,  dimensionless 

m  exponent  in  tlie  expression  T)  —  K  — —  y,  dimensionless 

N  exponent  in  the  expression  (y/s)1^,  dimensionless 
n  index  of  refraction,  dimensionless 

Pr  Prandtl  number,  dimensionless 

p  static  pressure,  slugs/ sq  ft 

R  universal  gas  constant,  ft/°R 

Upgx 

Re  local  Reynolds  number  — - ,  dimensionless 


r  temperature  ratio  — — ,  dimensionless 

S  interferogram  fringe  shift,  dimensionless 

s  distance  along  light  path,  ft 
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T  gas  static  temperature  at  point  x,  y,  z  in  the  "boundary  layer,  °R 
T*  gas  total  temperature,  °R 

U  velocity  of  gas  at  outer  edge  of  the  boundary  layer,  ft/sec 
u  velocity  of  gas  in  x  direction,  ft/sec 
v  velocity  of  light  at  point  x,  y,  z;  ft/sec 

x  distance  along  surface  measured  from  starting  point  of  boundary 
layer,  ft 

y  distance  normal  to  the  surface,  ft 

z  distance  along  span  of  test  section,  ft 

3  exponent  in  the  specific  heat-temperature  relation  cp/cp  5  =  (T/T§)^ 

\x  absolute  viscosity,  slugs/( sec)  (ft) 

I  |-(kp5)  K2Re2m  ( ?  dimensionless 

p  density,  slugs/cu  ft 

T  thickness  of  glass  vail,  ft 

T*  - -I_ — -  ft 

cos  (e2/ngr 

0)  exponent  in  viscosity  -  temperature  relation,  M-/ “  (T/Tg)*0, 
dimensionless 

Subscripts: 

1 

a  apparent  value 

f  inside  glass-wall  boundary-layer  film 

g  glass  side  wall 

l  laminar  boundary  layer 

o  compensating  chamber  conditions 

s  condition  at  border  of  laminar  sublayer  and  turbulent  region 

t  turbulent  boundary  layer 
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r 


c 


<• 


X 

n 

e 

0 

v 

w 

6 

G 

1 


°p 

ratio  of  specific  heats,  — ,  dimensionless 

°v 


thickness  of  thermal  boundary  layer,  ft 

ratio  of  distance  between  glass  and  object  plane  to  distance 
between  the  glass  and  test  section  centerline,  dimensionless 

ratio  of  distance  between  position  two,  figure  3,  and  object  plane 
to  distance  between  position  two  and  test  section  centerline, 
eL/ 2  -  6f 

y- - — ,  dimensionless 


dimensionless 


K  — —  y,  dimensionless 
x  J 

-^L*/z)  2. — - ,  dimensionless 

K  (kp&)  Rem  (L*/x)  d£ 


slope  of  ray  path  at  point  y,  z;  dimensionless 

wavelength  of  light  used,  ft 

vacuum 

wall 

condition  at  border  of  boundary  layer  and  free  stream 

referring  to  optical  path  from  inside  of  exit  window  to  the  image 
plane 

position  at  which  light  ray  enters  boundary  layer  under  investiga¬ 
tion 

position  at  which  light  ray  leaves  the  boundary  layer  under  invest! 
gation 


* 


* 


2 
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APPENDIX  B 


DERIVATION  OF  EQUATIONS  (l)  THROUGH  (8) 


The  equation 


y"  1  dv  .  1  dv  , 

1  +  yi2  “v5y  v3zy 


describes  the  trajectory  of  a  light  ray  in  a  nonhomogene ous  medium 
(reference  7).  By  the  following  substitutions,  equation  (Bl)  can  be  put 
in  a  more  useful  form: 


y*  =  tan  9  «  9 


n  =  kp  +  1 


(for  air) 


then 


kp  «  1 


dz  dy 


When  this  equation  is  integrated  with  the  houndary  condition  z  =  0 
when  0  =  0  yields 


e  =  k  z 
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If  the  equations 


> 


y 


are  substituted  into  equation  (B4), 

$  .  -  K*  (*ps)  $ 


dr 

dn 


By  integrating  equation  (B6)  between  the  limits  tj  —  tji  when  ( 
and  T|  =  tj  when  5=5  and  by  substituting 

2 


25  2  =  K2  (kps)  Re2*/^- 


) 


and 


« -  ^ 


then 


i  =  - 


dn 


dr/ dT) 


This  equation  represents  the  light-ray  trajectories  in  terms  of 
variables* 

By  letting  0  =  it  is  found  that 

d| 

K  (kp5)  Rem  (L’/x) 

6  =  - 


(B5) 


(B6) 
=  0 

(B7) 

(B8) 


(1) 

the  new 


(B9) 


0 
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If  the  object  plane  is  considered  to  be  at  a  distance 
L1 

€•  —  from  position  2  in  the  test  section 

u 


ya  y2  -  C'  2  &2 


(BIO) 


•where  9  is  a  negative  angle. 
Substitution  of  the  new  variables  yields 


02 

E  —  e* 
®2  e 


(Bii) 


and 


C'  ‘*82 


(2) 


The  relation  generally  used  between  index  of  refraction  change  and 
fringe  shift  (reference  l)  is 


-  nQ) 


(3) 


The  optical -path  difference  between  a  light  path  through  the  test  sec¬ 
tion  and  a  corresponding  light  path  through  the  compensating  section 
(reference  4)  is 


SX  =  2n„  6„  +  n  ds  +  n  T*  +  n  (h  +  T )  -  2n  5 „  - 

f  f  I  g  o'  of 


noL’  -  v  -no(^rs;) 


(B12) 


By  combining  equations  (3)  and  (B12) ,  by  substituting  for  the  index 
of  refraction 


n  =  kp  +  1 


(4) 
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and  by  solving  for  Lkp_  it  is  found  that 

& 


Lkpa  «  2kpf  Sf  +  k  /  p  ds  +  -  /  9“  dz  + 


rigA  +  /  dz  -  L 


(B13) 


where 


Lo  _  i  “o  (  h  +  T 


A  =  T«  -  T  +  —  h+T  -  —  - 5- 

ng  V  /  ng\^cos02 


(B14) 


and  where 


^  ds  =y^  (i  +  e2)  dz  (i  +  |e2) 


(B15) 


when  powers  of  0  equal  to  four  or  greater  are  neglected. 

If  equation  (B13)  is  solved  for  pa/p1  it  can  be  shown  that 


—  =  2  — i  —  /  p  ds  +  |  /  02  dz  +  ^  (B16) 

px  px  L  L  px  /  2  kLpL  /  kp1  L 


When  the  temperature  ratio  r  =  T/Tg  and  the  law  of  state 
p  =  p/RrTg  are  introduced,  it  follows  that 


Pa  rl  8f  1  /  dz  1_  _  7  0  2  dz  . 

P^-2^T+Lrl  J  T+2Lrl7  —  + 


1  1 


Lg  „  A 


~  rn  /  0“  dz  +  r*.  r  = 
2L  1  J  kpg  L 


(B17) 


The  expression  for  n  can  "be  put  in  the  following  useful  form: 

© 
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By  substituting  the  new  variables  into  equation  (B24),  it  can  be 
shcrwn  that 


Pa  ^  (jfjf 

1  +ri  2^7 


2  WO 


-aU 

r/f  Jo 


/T®  dt  “  2  rl  ^2  ®2 


(B25) 


where 


pa\  f  £a\  ,  f  £a\ 

plJf+\pl/l-2  VPl/« 


£a\  „  2  fl  ^ 

pJf  rf  L 


1  2,/T 


62  ft  d| 


« 

1/  “  "  2 


rl  1 2  ©2‘ 
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TABLE  I  -  RESULTS  OF  NUMERICAL  EXAMPLE 

[Assumptions:  tj^  determined  at  dr/drj  a  0.9(dr/dT])w  for  laminar 
boundary  layer ;  T]]_  »  0.01  for  turbulent  boundary  layer;  minimum 

y^  =  0.004  in.;  maximum  Re  =  5X10^;  pressure  =  26  in.  Hg  vacuum; 
temperature  a  70°  F;  length  of  cooled  wall  a  channel  length; 
k  a  0.1166  ft/slug;  r^  a  l7] 


1 

2 

3 

4 

5 

6 

7 

8 

9 

r 

w 

Mini¬ 
mum  X 
(ft) 

Maxi¬ 
mum  L 
(in.) 

Maximum 

velocity 

U 

(ft/sec) 

e  at 

ya 

1 

yl 

Pl-Pa 

pl 

^w-6 

x/L 

(Pa/Pl)f 

Pa/pl" 

Laminar 

0.343 

3.0 

1.7 

210 

0.78 

3.3 

0.13 

boundary 

.628 

1.2 

3.1 

540 

.79 

El 

2.5 

4.5 

.04 

layer 

.890 

1.0 

6.8 

630 

.86 

1.2 

1.8 

.02 
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Figure  1.  -  Schematic  diagram  ot  path  or  light  ray  through  two-dimensional  test  section. 
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(b)  Turbulent  "boundary  layer. 
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Figure  7.  -  Effect  of  wall  "boundary  layer  on 
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Concluded.  Effect  of  refraction 
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